The evolution of the dust temperatures of galaxies in the SFR–M_∗ plane up to z ~ 2

We study the evolution of the dust temperature of galaxies in the SFR− M∗ plane up to z ~ 2 using far-infrared and submillimetre observations from the Herschel Space Observatory taken as part of the PACS Evolutionary Probe (PEP) and Herschel Multi-tiered Extragalactic Survey (HerMES) guaranteed time key programmes. Starting from a sample of galaxies with reliable star-formation rates (SFRs), stellar masses (M_∗) and redshift estimates, we grid the SFR− M_∗parameter space in several redshift ranges and estimate the mean dust temperature (T_(dust)) of each SFR–M_∗ − z bin. Dust temperatures are inferred using the stacked far-infrared flux densities (100–500  μm) of our SFR–M_∗ − z bins. At all redshifts, the dust temperature of galaxies smoothly increases with rest-frame infrared luminosities (L_(IR)), specific SFRs (SSFR; i.e., SFR/M_∗), and distances with respect to the main sequence (MS) of the SFR− M_∗ plane (i.e., Δlog (SSFR)_(MS) = log [SSFR(galaxy)/SSFR_(MS)(M_∗,z)]). The T_(dust) − SSFR and T_(dust) − Δlog (SSFR)_(MS) correlations are statistically much more significant than the T_(dust) − LIR one. While the slopes of these three correlations are redshift-independent, their normalisations evolve smoothly from z = 0 and z ~ 2. We convert these results into a recipe to derive T_(dust) from SFR, M_∗ and z, valid out to z ~ 2 and for the stellar mass and SFR range covered by our stacking analysis. The existence of a strong T_(dust) − Δlog (SSFR)_(MS) correlation provides us with several pieces of information on the dust and gas content of galaxies. Firstly, the slope of the T_(dust) − Δlog (SSFR)_(MS) correlation can be explained by the increase in the star-formation efficiency (SFE; SFR/M_(gas)) with Δlog (SSFR)_(MS) as found locally by molecular gas studies. Secondly, at fixed Δlog (SSFR)_(MS), the constant dust temperature observed in galaxies probing wide ranges in SFR and M_∗ can be explained by an increase or decrease in the number of star-forming regions with comparable SFE enclosed in them. And thirdly, at high redshift, the normalisation towards hotter dust temperature of the T_(dust) − Δlog (SSFR)_(MS) correlation can be explained by the decrease in the metallicities of galaxies or by the increase in the SFE of MS galaxies. All these results support the hypothesis that the conditions prevailing in the star-forming regions of MS and far-above-MS galaxies are different. MS galaxies have star-forming regions with low SFEs and thus cold dust, while galaxies situated far above the MS seem to be in a starbursting phase characterised by star-forming regions with high SFEs and thus hot dust

The deepest Herschel-PACS far-infrared survey: number counts and infrared luminosity functions from combined PEP/GOODS-H observations

We present results from the deepest Herschel-Photodetector Array Camera and Spectrometer (PACS) far-infrared blank field extragalactic survey, obtained by combining observations of the Great Observatories Origins Deep Survey (GOODS) fields from the PACS Evolutionary Probe (PEP) and GOODS-Herschel key programmes. We describe data reduction and theconstruction of images and catalogues. In the deepest parts of the GOODS-S field, the catalogues reach 3σ depths of 0.9, 0.6 and 1.3 mJy at 70, 100 and 160 μm, respectively, and resolve ~75% of the cosmic infrared background at 100 μm and 160 μm into individually detected sources. We use these data to estimate the PACS confusion noise, to derive the PACS number counts down to unprecedented depths, and to determine the infrared luminosity function of galaxies down to L_(IR) = 10^(11) L⊙ at z ~ 1 and L_(IR) = 10^(12) L⊙ at z ~ 2, respectively. For the infrared luminosity function of galaxies, our deep Herschel far-infrared observations are fundamental because they provide more accurate infrared luminosity estimates than those previously obtained from mid-infrared observations. Maps and source catalogues (>3σ) are now publicly released. Combined with the large wealth of multi-wavelength data available for the GOODS fields, these data provide a powerful new tool for studying galaxy evolution over a broad range of redshifts

The PEP Survey: Infrared Properties of Radio-Selected AGN

By exploiting the VLA-COSMOS and the Herschel-PEP surveys, we investigate the Far Infrared (FIR) properties of radio-selected AGN. To this purpose, from VLA-COSMOS we considered the 1537, F[1.4 GHz]>0.06 mJy sources with a reliable redshift estimate, and sub-divided them into star-forming galaxies and AGN solely on the basis of their radio luminosity. The AGN sample is complete with respect to radio selection at all z<~3.5. 832 radio sources have a counterpart in the PEP catalogue. 175 are AGN. Their redshift distribution closely resembles that of the total radio-selected AGN population, and exhibits two marked peaks at z~0.9 and z~2.5. We find that the probability for a radio-selected AGN to be detected at FIR wavelengths is both a function of radio power and redshift, whereby powerful sources are more likely to be FIR emitters at earlier epochs. This is due to two distinct effects: 1) at all radio luminosities, FIR activity monotonically increases with look-back time and 2) radio activity of AGN origin is increasingly less effective at inhibiting FIR emission. Radio-selected AGN with FIR emission are preferentially located in galaxies which are smaller than those hosting FIR-inactive sources. Furthermore, at all z<~2, there seems to be a preferential (stellar) mass scale M ~[10^{10}-10^{11}] Msun which maximizes the chances for FIR emission. We find such FIR (and MIR) emission to be due to processes indistinguishable from those which power star-forming galaxies. It follows that radio emission in at least 35% of the entire AGN population is the sum of two contributions: AGN accretion and star-forming processes within the host galaxy.Comment: 13 pages, 14 figures, to appear in MNRA

Molecular Gas, Dust and Star Formation in Galaxies: II. Dust properties and scalings in \sim\ 1600 nearby galaxies

We aim to characterize the relationship between dust properties. We also aim to provide equations to estimate accurate dust properties from limited observational datasets. We assemble a sample of 1,630 nearby (z<0.1) galaxies-over a large range of Mstar, SFR - with multi-wavelength observations available from wise, iras, planck and/or SCUBA. The characterization of dust emission comes from SED fitting using Draine & Li dust models, which we parametrize using two components (warm and cold ). The subsample of these galaxies with global measurements of CO and/or HI are used to explore the molecular and/or atomic gas content of the galaxies. The total Lir, Mdust and dust temperature of the cold component (Tc) form a plane that we refer to as the dust plane. A galaxy's sSFR drives its position on the dust plane: starburst galaxies show higher Lir, Mdust and Tc compared to Main Sequence and passive galaxies. Starburst galaxies also show higher specific Mdust (Mdust/Mstar) and specific Mgas (Mgas/Mstar). The Mdust is more closely correlated with the total Mgas (atomic plus molecular) than with the individual components. Our multi wavelength data allows us to define several equations to estimate Lir, Mdust and Tc from one or two monochromatic luminosities in the infrared and/or sub-millimeter. We estimate the dust mass and infrared luminosity from a single monochromatic luminosity within the R-J tail of the dust emission, with errors of 0.12 and 0.20dex, respectively. These errors are reduced to 0.05 and 0.10 dex, respectively, if the Tc is used. The Mdust is correlated with the total Mism (Mism \propto Mdust^0.7). For galaxies with Mstar 8.5<log(Mstar/Msun) < 11.9, the conversion factor \alpha_850mum shows a large scatter (rms=0.29dex). The SF mode of a galaxy shows a correlation with both the Mgass and Mdust: high Mdust/Mstar galaxies are gas-rich and show the highest SFRs.Comment: 24 pages, 28 figures, 6 tables, Accepted for publication in A&

IRAC Excess in Distant Star-Forming Galaxies: Tentative Evidence for the 3.3μ\mum Polycyclic Aromatic Hydrocarbon Feature ?

We present evidence for the existence of an IRAC excess in the spectral energy distribution (SED) of 5 galaxies at 0.6<z<0.9 and 1 galaxy at z=1.7. These 6 galaxies, located in the Great Observatories Origins Deep Survey field (GOODS-N), are star forming since they present strong 6.2, 7.7, and 11.3 um polycyclic aromatic hydrocarbon (PAH) lines in their Spitzer IRS mid-infrared spectra. We use a library of templates computed with PEGASE.2 to fit their multiwavelength photometry and derive their stellar continuum. Subtraction of the stellar continuum enables us to detect in 5 galaxies a significant excess in the IRAC band pass where the 3.3 um PAH is expected. We then assess if the physical origin of the IRAC excess is due to an obscured active galactic nucleus (AGN) or warm dust emission. For one galaxy evidence of an obscured AGN is found, while the remaining four do not exhibit any significant AGN activity. Possible contamination by warm dust continuum of unknown origin as found in the Galactic diffuse emission is discussed. The properties of such a continuum would have to be different from the local Universe to explain the measured IRAC excess, but we cannot definitively rule out this possibility until its origin is understood. Assuming that the IRAC excess is dominated by the 3.3 um PAH feature, we find good agreement with the observed 11.3 um PAH line flux arising from the same C-H bending and stretching modes, consistent with model expectations. Finally, the IRAC excess appears to be correlated with the star-formation rate in the galaxies. Hence it could provide a powerful diagnostic for measuring dusty star formation in z>3 galaxies once the mid-infrared spectroscopic capabilities of the James Webb Space Telescope become available.Comment: 25 pages, 4 figures, accepted by Ap

Total Molecular Gas Masses of Planck - Herschel Selected Strongly Lensed Hyper Luminous Infrared Galaxies

We report the detection of CO(1 - 0) line emission from seven Planck and Herschel selected hyper luminous (LIR(8-1000um) > 10^13Lsun) infrared galaxies with the Green Bank Telescope (GBT). CO(1 - 0) measurements are a vital tool to trace the bulk molecular gas mass across all redshifts. Our results place tight constraints on the total gas content of these most apparently luminous high-z star-forming galaxies (apparent IR luminosities of LIR > 10^(13-14) Lsun), while we confirm their predetermined redshifts measured using the Large Millimeter Telescope, LMT (zCO = 1.33 - 3.26). The CO(1 - 0) lines show similar profiles as compared to Jup = 2 -4 transitions previously observed with the LMT. We report enhanced infrared to CO line luminosity ratios of = 110 (pm 22) Lsun(K km s^-1 pc^-2)^-1 compared to normal star-forming galaxies, yet similar to those of well-studied IR-luminous galaxies at high-z. We find average brightness temperature ratios of = 0.93 (2 sources), = 0.34 (5 sources), and = 0.18 (1 source). The r31 and r41 values are roughly half the average values for SMGs. We estimate the total gas mass content as uMH2 = (0.9 - 27.2) x 10^11(alphaCO/0.8)Msun, where u is the magnification factor and alphaCO is the CO line luminosity to molecular hydrogen gas mass conversion factor. The rapid gas depletion times are, on average, tau = 80 Myr, which reveal vigorous starburst activity, and contrast the Gyr depletion timescales observed in local, normal star-forming galaxies.Comment: published in MNRAS, 9pages, 5fig

Automated mining of the ALMA archive in the COSMOS field (A3COSMOS): II. Cold molecular gas evolution out to Redshift 6

We present new measurements of the cosmic cold molecular gas evolution out to redshift 6 based on systematic mining of the ALMA public archive in the COSMOS deep field (A3COSMOS). Our A3COSMOS dataset contains ~700 galaxies (0.3 < z < 6) with high-confidence ALMA detections in the (sub-)millimeter continuum and multi-wavelength spectral energy distributions (SEDs). Multiple gas mass calibration methods are compared and biases in band conversions (from observed ALMA wavelength to rest-frame Rayleigh-Jeans(RJ)-tail continuum) have been tested. Combining our A3COSMOS sample with ~1,000 CO-observed galaxies at 0 < z < 4 (75% at z < 0.1), we parameterize galaxies' molecular gas depletion time and molecular gas to stellar mass ratio (gas fraction) each as a function of the stellar mass, offset from the star-forming main sequence (Delta MS) and cosmic age (or redshift). Our proposed functional form provides a statistically better fit to current data (than functional forms in the literature), and implies a "downsizing" effect (i.e., more-massive galaxies evolve earlier than less-massive ones) and "mass-quenching" (gas consumption slows down with cosmic time for massive galaxies but speeds up for low-mass ones). Adopting galaxy stellar mass functions and applying our function for gas mass calculation, we for the first time infer the cosmic cold molecular gas density evolution out to redshift 6 and find agreement with CO blind surveys as well as semi-analytic modeling. These together provide a coherent picture of cold molecular gas, SFR and stellar mass evolution in galaxies across cosmic time

The Far-Infrared, UV and Molecular Gas Relation in Galaxies up to z=2.5

We use the infrared excess (IRX) FIR/UV luminosity ratio to study the relation between the effective UV attenuation (A_IRX) and the UV spectral slope (beta) in a sample of 450 1<z<2.5 galaxies. The FIR data is from very deep Herschel observations in the GOODS fields that allow us to detect galaxies with SFRs typical of galaxies with log(M)>9.3. Thus, we are able to study galaxies on and even below the main SFR-stellar mass relation (main sequence). We find that main sequence galaxies form a tight sequence in the IRX--beta plane, which has a flatter slope than commonly used relations. This slope favors a SMC-like UV extinction curve, though the interpretation is model dependent. The scatter in the IRX-beta plane, correlates with the position of the galaxies in the SFR-M plane. Using a smaller sample of galaxies with CO gas masses, we study the relation between the UV attenuation and the molecular gas content. We find a very tight relation between the scatter in the IRX-beta plane and the specific attenuation (S_A), a quantity that represents the attenuation contributed by the molecular gas mass per young star. S_A is sensitive to both the geometrical arrangement of stars and dust, and to the compactness of the star forming regions. We use this empirical relation to derive a method for estimating molecular gas masses using only widely available integrated rest-frame UV and FIR photometry. The method produces gas masses with an accuracy between 0.12-0.16 dex in samples of normal galaxies between z~0 and z~1.5. Major mergers and sub-millimeter galaxies follow a different S_A relation.Comment: 11 pages, 6 pages appendix, 11 figures, accepted to Ap

Probing the Timescale of the 1.4 GHz Radio emissions as a Star formation tracer

Radio used as a star formation rate (SFR) tracer presents enormous advantages by being unaffected by dust and radio sources being pinpointed at the sub-arc-second level. The interpretation of the low frequency 1.4 GHz luminosity is hampered by the difficulty in modeling the cosmic ray paths in the interstellar medium, and their interactions with the magnetic field. In this work, we compare the SFR derived from radio observations, and the ones derived from spectral energy distribution (SED) modeling. We aim at better understand the behavior of the SFR radio tracer, with a specific emphasis on the link with star-formation histories. We used the SED modeling code Code Investigating GALaxy Emission, CIGALE, with a non-parametric star formation history model (SFH) and fit the data over the wavelength range from the ultraviolet (UV) up to the mid-infrared (mid-IR). We interpret the difference between radio and SED-based SFR tracers in the light of recent gradients in the derived SFH. To validate the robustness of the results, we checked for any remaining active galaxy nuclei (AGN) contribution and tested the impact of our SFH modeling approach. Approximately 27% our galaxies present a radio SFR (SFR$_{\rm radio}$) at least ten times larger than the instantaneous SFR from SED-fitting (SFR$_{\rm SED}$). This trend affects primarily the galaxies that show a declining SFH activity over the last 300 Myr. Both SFR indicators converge toward a consistent value, when the SFHs are averaged over a period larger than 150 Myr to derive SFR$_{\rm SED}$. Although the radio at low frequency 1.4 GHz is a good tracer of the star formation activity of galaxies with constant or increasing SFH, our results indicate that this is not the case for galaxies that are quenching. Our analysis suggests that the star formation time sensitivity of the radio low frequency could be longer than 150 Myr.Comment: 10 pages, 10 figure